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Abstract. We calculate the average of two characteristic polynomials for the 
real Ginibre ensemble of asymmetric random matrices, and its chiral counterpart. 
Considered as quadratic forms they determine a skew-symmetric kernel from which 
all complex eigenvalue correlations can be derived. Our results are obtained in a very 
simple fashion without going to an eigenvalue representation, and are completely new 
' in the chiral case. They hold for Gaussian ensembles which are partly symmetric, with 

kernels given in terms of Hermite and Laguerre polynomials respectively, depending 
on an asymmetry parameter. This allows us to interpolate between the maximally 
CS| ' asymmetric real Ginibre and the Gaussian Orthogonal Ensemble, as well as their chiral 

\ counterparts. 

PACS numbers: 02.10.Yn, 0250.-r, 0540.-a 
95 ■ 1- Introduction 

o 

Random Matrix Theory is known to enjoy a wide range of applications in the physical 
^ \ sciences and beyond. This remains true when the eigenvalues of the operator to 
\ be described move into the complex plane. However, the ensemble that is perhaps 
the most interesting of these, the real Ginibre ensemble p] dealing with real-valued 
asymmetric matrix entries, has turned out to be the most difficult. Possible applications 
of these ensembles include neural networks [2], directed Quantum Chaos [3], Quantum 
Chromodynamics [1], financial markets [5], and quantum information theory [6]. 

The mathematical difficulty in solving these ensembles is due to the fact that they 
allow for combinations of both real and complex conjugate eigenvalue pairs, with their 
characteristic equation having only real entries. Apart from results on the spectral 
density [Zl El [3] an eigenvalue representation [HI HO] was derived as a starting point for 
studying systematically higher order eigenvalue correlation functions. Only very recently 
was their Pfaffian structure explicitly revealed [HI [12], and the probability p7v,A; that an 
N X N matrix has exactly k real eigenvalues [HI [13] as well as all correlations for k = 
were computed [13j. Finally the complete solution for all real and complex eigenvalue 
correlations was achieved independently by three different groups [TU [IS, [161 [13 [IE] • 
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In this paper we will give a very simple derivation for the generating kernel of 
all complex eigenvalue correlations in the general so called elliptic case, dealing with 
partly symmetric matrices depending on an asymmetry parameter. We also present 
new results for the chiral real Ginibre ensemble as a two-matrix model which has not 
yet been considered. 

In the next section we will explain the relation between the complex eigenvalue 
density and characteristic polynomials. After defining the elliptic real Ginibre ensemble 
and its new chiral extension we consider their corresponding results in two separate 
sections |3] and HI Our conclusions are presented in section [51 



2. The role of characteristic polynomials generating the kernel 

We start with the simplest ensemble considered here, the real Ginibre ensemble at 
maximal asymmetry. It is just given by a Gaussian measure in the space of real, 
asymmetric N x N matrices which is invariant under orthogonal transformations: 

- n (^) - DA e-^-^^^ (1) 

i,j = l ^ * / \ / 

The eigenvalues Aj obey the equation det[Aj — A] =0, and thus are real or occur in 
complex conjugate pairs. They enjoy the following ordered joint probability density 
function (jpdf ) [10] 

N N 

d^{X,, A2, . . . , A^) =CM-dX,... dXM ■ 1[{X, - A,) • J] f{Xk) , (2) 

i<j k 

with some positive definite weight function /(A^) = f(Xk) and a normalisation constant 
Cat. Here the eigenvalues are ordered as follows if they are real: Ai > A2 > . . . , if 
they are complex: ReAi = ReA2 > ReAs = ReA4 > . . . , ImAi = — ImA2 > 0, ImAs = 
— ImA4 > 0, . . . , and similarly if they are mixed (see also [T8]). 

This implies that the spectral density of complex eigenvalues of the + 2 
dimensional ensemble, which can be obtained by inserting a two-dimensional delta- 
function in the complex plane, is proportional to 

R%+2,iW i{X - A)/(A)2(det[A - A] det[A - A])n • (3) 

The brackets mean the average over the ensemble ([1]) with the partition function 
Zn = jBAe-'2^^^^. There is an additional contribution to the total spectral density 
from the real eigenvalues, which is obtained by inserting a delta-function on the real 
axis, which we do not consider here. We are therefore led to consider the following 
correlation of two characteristic polynomials of A 

^^(A,7) = (det[A - A] det[7 -A^])n = ^\ with A^7. (4) 

A - 7 

\ The Hermitian conjugate is put here merely to stress the analogy to other ensembles discussed below, 
as for real matrices det[A'^] = det[A]. 
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It determines the antisymmetric kernel /C^~^(A,7), from which all correlation functions 
of complex eigenvalues follow for even A^. While the result for odd is known in 
the real Ginibre ensemble [18] , very recently a general technique has been proposed for 
obtaining the odd result from even by removing an eigenvalue [T9|. The kernel in 
eq. @ has been derived in [H] using the Edelman result [lO] for the density of complex 
eigenvalues, and using eq.([3]). Here we will give an independent derivation which makes 
clear why the result is so simple, when the jpdf eq. ([2]) is so complicated. 

The relation eq. (jl]) is far more general. Not only does it hold for the other real 
Ginibre ensembles to be introduced below, but it also holds for other symmetry classes 
with complex eigenvalues having unitary or symplectic invariance. For the quaternionic 
Ginibre ensembles at /5 = 4 an identical relation to eq. (jl]) was shown in [20j to give 
the skew-symmetric kernel. In the Ginibre ensembles with unitary symmetry (3 = 2 the 
kernel is symmetric and the following modified, simpler relation is known to hold [21] 

(det[A - A] det[7 - ^t])^ = /C^='(A,7) • (5) 

The argument we just presented above for the real Ginibre ensemble at maximal 
asymmetry can easily be translated to the partially symmetric case depending on an 
asymmetry parameter r. Here in the large- A^ limit the complex eigenvalues lie inside 
an ellipse with axes ~ (1 ± r) [7]. It is known [9] that its jpdf is related to eq. ([2]) by a 
simple rescaling of the eigenvalues, and we readily obtain 

J-^(A,7;r) = (det[A - {S + vA)] det[7 - {S + vAf])j, ^ '}^- ^' , (6) 

A — 7 



with 



rG[0,l], v' = \-^. (7) 
1 + r 



The average is with respect to the following partition function 
Zn = [bS da exp 



^ TiiSS^ + AA^ 



(8) 



2(1 + r 

and we consider the eigenvalues of the partly symmetric matrix J = S + vA. Here 
S and A are N x N matrices being symmetric and antisymmetric respectively, with a 
particular choice of variance. The limiting case r = brings us back to the ensemble 
eq. ([I]) while setting r = 1 would lead to the Gaussian Orthogonal Ensemble. However, 
in that case the eigenvalues become real and this limit is subtle. 

The second ensemble we consider in this paper is the chiral counterpart of the real 
Ginibre ensemble. Following [1] and its extension to a two-matrix model [22] we define 
the following chiral real Ginibre ensemble {ch) with a particular variance n, 

I r 71 
Z'f!; = DADB exp Ti{AA^ + BB^) . (9) 



2 

Again we compute the average of characteristic polynomials to obtain the kernel 



^^'^(A,7;yu) = (det[A-M]det[7-M^])^ = ^ ^ ' , (10) 



A2 -72 

M = I " „^ " '^"^^ 1 . fill 
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Here both A and B are rectangular x (A^ + u) matrices without further symmetry 
among the real matrix elements. They are drawn independently from the ensemble 
([1]) extended to u > 0. The asymmetry parameter is given here by n G [0,1], where 
^ = 1 denotes maximal asymmetry, and /j, = takes us back to the chiral Gaussian 
Orthogonal Ensemble. In [4J initially a one-matrix model was proposed, replacing B 
by the identity. Whilst we expect that in the large-A^ limit both lead to the same 
universal result our choice allows for an eigenvalues basis, as in the corresponding chiral 
extensions of Ginibre at /3 = 2 [22] and (3 = 4 [23], having complex and quaternion 
real matrix elements, respectively. The ensemble in eq. (Q has not been solved before 
and we will give a completely new result below, depending parametrically on u. For 
maximal asymmetry it corresponds to class 2P in [23] for real elements. 
The 2N + u eigenvalues Aj of the matrix M defined in eq. ( fTTI) satisfy 

= det[A-M] = X^detlX^ -{A + fiB){A^ - fiB^)] . (12) 

It will be shown elsewhere that the jpdf of these eigenvalues is again of the form in eq. 
([2]). From that it follows that the kernel derived from JF^(A,7; fi) again determines all 
correlation functions of complex eigenvalues. 

A peculiarity of the chiral ensemble is the following: the non-zero eigenvalues 
solving the second equation in ( |T2|) are real but not necessarily positive. Thus the 
eigenvalues of M can have both real and purely imaginary eigenvalues as well as complex 
conjugate eigenvalue pairs. Moreover, all non-zero eigenvalues come in ± pairs due to 
the chirality of the matrix M. 



3. Characteristic polynomials for the real Ginibre ensemble 

For pedagogical reasons we begin with the maximally asymmetric case eq. (jlj). The 
partly symmetric case at r 7^ is given as a second example below. 

^^(A, 7) = ^ / DA e-^T--^^^ det[A - A] det[7 - A^] . (13) 
Zn J 

Writing the determinants in terms of two A^- dimensional complex Grassmann vectors 
r]i and Q, with z = 1, . . . , A^, we obtain 



^7v(A,7) 



1 
2n 



DA / drj exp 



dQ drj exp 



(14) 



after integrating out the Gaussian matrix A. Here and in the following we will use 
summation conventions over double indices. The last term in the exponent can be 
written as 



(15) 
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With the help of a complex Hubbard- Stratonovich (HS) transformation we can 
bilinearise and integrate out the Grassmann variables: 



^7v(A,7) 



^ j j dC dri exp [-\z\'^ - XQQ - -fri*r]i + zQt]* + zQr]^] 



N 



d'z e-l^l' (At+I^H 



2\N 



(A7)" 



n=0 



(16) 



This gives a polynomial with leading power (A7)^ as expected. Thus our first main 
result leads to the following antisymmetric kernel 



N 



/Cjv(A,7) = {X-i)N\J2 



n=0 



n! 



;i7) 



which is enough to derive all complex correlation functions. On setting 7 = A and 
multiplying by the weight /(A)^, Edelman's complex density [8] in terms of an incomplete 
exponential follows. It is remarkable that it only depends on | while the jpdf eq. ([2]) is 
not isotropic. We note that Edelman derived his result using methods from multivariate 
statistics, and not from the jpdf. 

We now turn to the partly symmetric case with r G [0, 1], where we can follow the 
same path, 



•^iv(A,7;r) 



■■^ [dSDA e-^^(^^''+^^'') det[A - (5 + vA)] det[7 - {S + vAf 
Zn J 



1 



D^DA / dC dr] exp 



-(4-A^)-AC0-7^N. 



2(1 + r 

+ Ci{Sij + vAij)Q + inil{Sij - vAij)r]j 



After symmetrising and antisymmetrising the terms in the last line, e.g. QSijCj = 
\Sij{QC,j + CjO); "we can complete the squares in Sij and Aij respectively, and integrate 
them out to obtain 



J'7v(A,7;r) 



dQ drj exp 



ACO - IVH - cl{{Qm*Q + {riHMm)) 



Here we have introduced the constants c\ = |(1 + r)(l ± t>^). The quartic terms in the 
Grassmann variables can be rewritten using two real HS transformations for the first 
line of eq. f[T^ . and two complex ones for the second line: 
1 



-^7v(A,7;^) 



TX--' 



dx dy / d z d w / dC, dr] exp — x — y — C*(A + 2ic-x)Q 



r]*{-f + 2ic-y)rij - 2\z\'^ - 2\w\^ + 2c+{zrijQ - zri*Q) - 2c^{wr]jC + wr]*Q) 



— dx dy d z d w exp 



-x-y- 2\z\ 



(A + 2ic_x)(7 + 2ic-y) + 4:cl\z\'^ + 4c^|if;| 



2\w\ 



N 



(20) 
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Expanding the last factor twice into binomial series in powers of \z\'^ and |wp we can 
apply the following integral representation of the Hermite polynomials: 



l—rt U—n \ V ^ ' / \ 



where using eq. ([7]) we have 2c_ = v2r. This eliminates the two real integrations. 
After integrating out the two remaining complex variables z and w we finally arrive at 

N I 

I / A \ / \ 

(22) 

l=Q k=0 

As a check this is again a polynomial with leading order (A7)^. Although our result eq. 
( |22|) could be further simplified this form is most useful for obtaining the antisymmetric 
kernel by applying the Christoffel-Darboux formula to the inner sum: 

ICU\r,r) = wfi(I) '"'(//,« (^)//,(-i=)-h« A)) . (23) 



1=0 

This coincides precisely with the kernel of skew-orthogonal Hermite polynomials derived 
in [T7] via the jpdf, which is much more elaborate. As was shown there independently, 
this kernel generates all complex eigenvalue correlation functions of the partly symmetric 
ensemble eq. ([H]) for even A^, depending parametrically on r. 

A similar kernel given in terms of orthogonal (for (3 = 2) [25] or skew-orthogonal (for 
(3 = 4) [26] Hermite polynomials is known for the partly symmetric Ginibre ensembles. 

4. Characteristic polynomials for the chiral real Ginibre ensemble 

In this section we present the calculation only for the partly symmetric case of the 
chiral extension of the real Ginibre ensemble, depending on asymmetry parameter /x. 
The simpler result at maximal asymmetry with /i = 1 is given at the end of this section. 

A -{A + ^B) 

-{A^ - fiB^) A 



m\r,f^) = ^ hADB e-t^(^-^^+^^^)det 



X det 



7 -i^-f^B) 
-{A^ + fiB^) 7 



■y'ch 



JbABB JdijdijdCd^ exp [ - ^{Al + Bl) - X{r]*r], + ra^a) - l{QQ + 



+ r^;(A, + ^^B,,)iJ, + ra{K^ - i^BDn, + q{a,, - /.s,,)^, + ^:(a^, + /.S^JOJ -(24) 

Here we have written each determinant of size 2A^+z/ in terms of two Grassmann vectors, 
rii[C,i) and ipaifa) of size and N + u, respectively. Our summation conventions imply 
for i = 1, . . . ,N, and for a = 1, . . . , N + u. In addition to fi we have a parameter n for 
the variance. After completing the square and integrating out the matrices A and B we 
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^NiK r,f^) = j drjdipdCdLp exp - X{r]*rii + tp^tpa) - 7(CCi + vl^a) 

Multiplying out and collecting all nonzero terms we need 6 new complex integrations 
to perform the HS transformations that bilinearise the Grassmann variables. We only 
give the result obtained after performing all Grassmann integrations: 

^^'^(A, 7; /i) =^J d\ d\ d^w d^p d'q d'z e-l«l^-l-I^H-l^-H^H5l^-l^l^ (26) 

X ((A — i6-u){'-y — i6-v) — + S'^pq)'^ ((A — i5-u){'j — i5-v) — + S'^pq)^'^'', 
where we have used the following abbreviations 

6l^-a±f,'). (27) 
n 



m - siwz + 5^M)^=f:H(-^>^)^-T(I)' 



Expanding the first factor in the second line of eq. fl26|) as 

with A7 = {X — i6-u){'-f — i6-v) , and likewise the second factor, we can use the following 
orthogonality relation 

- [ d^p e-IPl Vp' = 6kik\ . (29) 



Applying this first to the integrations over variables p and q, and then to w and z we can 
reduce the four sums to two. As a final step we employ the following complex integral 
representation for Laguerre polynomials 

i J d2«e-l"l'(A + m)'=(A + m)'=+'^ = k\ (-)'=A"L^(A') . (30) 

Whilst we did not find this representation in tables it can be easily verified from the 
standard representation of generalised Laguerre polynomials 

L^(x) - V (-^)"(^ + "^' X- (31) 



m=0 



Using the integral representation eq. (j30l) as well as its complex conjugate we finally 
arrive at the following result: 



1=0 ^ +^ k=0 



Ki\r,f^) = iV!(Ar + .)!5f(A7)^>:(^) } ] jjr^L^ ] LU ] ■ (32) 



It is a polynomial in A and 7 with the correct leading power (A7)^^"^'^. Looking back to 
the definition of the antisymmetric kernel in our chiral case, eq. ( |TOl) . we can read off the 
following result, after using the Christoffel-Darboux formula for Laguerre polynomials: 

N / e \ 4/ 



1=0 
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This gives our new kernel of the chiral real Ginibre ensemble, from which all its complex 
eigenvalue correlations follow. In particular for 7 = A it is proportional to a new complex 
eigenvalue density as in eq. ([3]). It is similar to the corresponding expressions for the 
kernel aX (3 = 2 [22] and /? = 4 [23j, also given in terms of Laguerre polynomials in the 
complex plane. 

After dealing with the general case we can go to maximal asymmetry, by setting 
fi = 1. In this limit only the leading power of the Laguerre polynomials contributes, 
and we obtain 

^-(A,,;,= l) = iV.(iV + (A,rg^^^(^) (34) 

for the characteristic polynomials, with lim^^i S"^ = ^. For the corresponding kernel we 
have to properly rescale with 5^ and we obtain 

/r)\2N N /^^2^ 

M.SlKr^Kr.,}-N<(N^^)<\l) (A7r(A^-7^)Efl^(^) .(35) 

When setting 7 = A and comparing to eq. ( ITTjl we again find a dependence on the 
modulus only, despite the anisotropic jpdf. Here the incomplete exponential is replaced 
by an incomplete modified /-Bessel function of the first kind. 

5. Conclusions 

We have calculated the expectation value of the product of two characteristic 
polynomials with respect the following two Gaussian random matrix models: the partly 
symmetric real Ginibre ensemble, and its chiral counterpart, a newly introduced two- 
matrix model. In our calculation we have used the supersymmetric method, without 
the need to explicitly go to an eigenvalue basis. In this simple way we can determine 
a skew-symmetric kernel which is the main building block for all complex eigenvalue 
correlation functions that can be written as Pfaffians. One could calculate this kernel 
directly from the joint eigenvalue distribution (jpdf), but this turns out to be a very 
difficult task. 

This kernel is given by a sum over Hermite polynomials for the real Ginibre case, 
depending on the asymmetry parameter. Here we have recovered a known, very recent 
result. In the chiral real Ginibre ensemble we find a new kernel given in terms of 
generalised Laguerre polynomials. In addition to the asymmetry ^ it depends on the 
parameter u labelling the number of exact zero eigenvalues. Our method offers an 
explanation of why the spectral density of complex eigenvalues is so simple, i.e. being 
an incomplete exponential or /-Bessel function at maximal asymmetry, while the jpdf 
is so complicated. 

One possible application of our new chiral result would be in field theory for Dirac 
operators with a real representation. The reason complex eigenvalues appear here is 
due to a chemical potential of the quarks. 
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It is an open question for the chiral ensemble if for all the kernel also determines 
the weight function /(A), and if both ingredients (i.e. kernel and weight) determine all 
correlation functions of real, complex and mixed eigenvalues. For the real Ginibre ensem- 
ble this fact is known to hold, and the similarity in structure makes this very suggestive. 
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